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Abstract 

Elastic proton scattering on 3 H nucleus is studied inbetween p- 3 H and n- 3 He thresholds, in 
the energy region where a-particles first excited state is inbeded in the continuum. For this aim 
Faddeev-Yakubovski equations are solved in configuration space, fully considering effects due to 
the isospin breaking as well as rigorously treating Coulomb interaction. Different realistic nuclear 
Hamiltonians are tested, elucidating open problems in nuclear interaction description. 
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I. INTRODUCTION 



During the last decade remarkable progress has been achieved in computational few- 
nucleon physics. First of all Green Function Monte-Carlo [l| and No Core Shell model 
techniques have been developed and perfectioned, enabling to describe bound nuclei contain- 
ing up to several nucleons. Lately these techniques have been also adopted to treat the sim- 



plest case of the elastic nucleon-nucleus scattering 



4|. Great progress have been achieved 



also in treating Coulomb interaction for three and four nucleon scattering , opening 
the new frontiers to test the nuclear interaction models and three-nucleon force (3NF) in 
particular. 

Regardless the continuous effort put to describe the strong nuclear interaction, it still 
remains the central issue in nuclear physics. If the on shell part of the nucleon-nucleon (NN) 
interaction is well constrained by the two nucleon scattering data, determination of the off 
shell properties and three-nucleon interaction in particular, which can be tested only in A>3 
nuclei, is highly untrivial task. As a consequence much effort have been devoted to describe 
the three nucleon system with special emphasis to N-d scattering problem. However due 
to the large size of the deuteron, it turns to be difficult to find observables, which reveal 
strong sensitivity to the off shell structure of NN interaction. Low energy nucleon-deuteron 
scattering - maybe apart the J 7r =l/2 + state, where strong correlations with triton binding 
energy are observed - is well described by NN interaction alone and is quite insensitive 
to its off shell structure. Only at large energies, well above deuteron break-up threshold, 
off-shell effects become more pronounced [?]. Unfortunately at large energies one faces 
increasing difficulty in controlling the relativistic effects and the large number of partial 
waves involved in the scattering process [§]] . Finding relatively simple system at low energy, 
which could be treated exactly and which is sensitive to the off-shell behavior of nuclear 
interaction is of great interest. Four nucleon continuum containing series of thresholds and 
resonances (see Fig. [1]) seems to present an ideal laboratory. Already the simplest case of 
n+ 3 H scattering, which is free of complications due to the Coulomb interaction, contains 
several broad resonances in the continuum and presents a serious test for nuclear interaction 
models: most of the realistic Hamiltonians by ~ 10% underpredict elastic cross section at 



the resonance region 



The most interesting but also the most complex four nucleon structure is continuum 



2 



of the 4 He, which contains numerous resonances and thresholds (see Fig. [T]). It has been 
suggested by Hofmann and Hale Q that 4 He system can be used database to fine 
tune the three-nucleon interaction. In particular delicate is the region between the p- 3 H 
and n- 3 He thresholds, due to the existence of the J"* = + resonance in the p- 3 H threshold 
vicinity. Accurate treatment of the charge symmetry breaking effects is required to separate 
the two thresholds: as example if one neglects the Coulomb repulsion between the two 



141 ]. becoming a 



protons the resonance of the a particle moves below the p- 3 H threshold 
bound state and changing completely the nearthreshold scattering dynamics. These facts 
reveal the particular relevance of this region in order to understand the isospin structure of 
the two and three nucleon interactions. 
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FIG. 1: The experimental spectra of ^iV bound and resonant states are presented, as obtained 
using R-matrix analysis [13]. In this figure resonances are indicated by the dashed lines. In the 
right pane of the figure the structure of 4 He is elucidated. 

This work aims to attract the attention to the p- 3 H elastic scattering problem, which is 
studied for the proton laboratory energies up to 1 MeV, i.e. below n- 3 He threshold. Six 
principally different realistic nuclear Hamiltonians will be used to study this system. 
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II. FADDEEV-YAKUBOVSKI EQUATIONS 



Four-nucleon problem is solved using Faddeev-Yakubovski (FY) equations in configura- 



tion space 



1 61 ] . In FY formalism four-particle wave function is written as a sum of 



18 amplitudes. From those we distinguish amplitudes of type K, which incorporate 3+1 
particle channels, while amplitudes of type H contain asymptotes of 2+2 particle channels. 
By interchanging order of the particles one can construct twelve different amplitudes of the 
type K and six amplitudes of the type H. 

Further we use the isospin formalism, i.e. we consider protons and neutrons as being 
degenerate states of the same particle - nucleon, having the mass fixed to h 2 jm = 41.47 
MeV-fm 2 . For the system of the identical particles one has only two independent FYA, 
one of type K and the other of type H . The other 16 FYA can be obtained by applying 
particle permutation operators {i.e. interchanging the order of the particles in the system). 



Similarly only two independent FY equations exist; by sin 



and including three-nucleon force, FY equations read: 17 



ing out K = K\ 2 3 and H = H 



31 
12 



18|: 



E - Ho - V n - JZ Vg ) K = V 12 (P + + P~) [(1 + Q)K + H} + l\V* + 




E-H -V 12 -y Vg \H = V 12 P[(1 + Q)K + H}. (1) 



i<j 



Here Vi 2 is the strong part of the NN interaction between nucleons (12), Vg is the Coulomb 
interaction between the nucleons % and j, the three-nucleon force (3NF) is represented by 
the terms V 2 \ and V^, which are symmetric for the cyclic particle permutations and which 
contain a part of 3NF acting in the particle cluster (123): Vi23=Vj| + V 2 \ + V^. The particle 
permutation operators P + , P~, P and Q we use are simply: 

p + = (p-y 1 = P23P12, 

Q = --^34, 

P — P13P24 — -P24-Pl3- 

Employing operators defined above, the total wave function of the four-nucleon is given by 

* = [1 + (1 + P + + P-)Q] (1 + P + + P-)K + (1 + P + + P")(l + P)H. (2) 

The numerical implementation of these equations is described in detail in 17| . This 
formalism enables to include Coulomb interaction as well as to test different realistic nuclear 
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interaction models, comprising non-local ones and models in conjunction with three-nucleon 
interaction. 

The 4 He system is predominantly total isospin T=0 state. However, the basis limited to 
T=0 does not allow unambiguous separation of the p- 3 H and n- 3 He channels. In order to 
separate these mirror channels we allow for full isospin-symmetry breaking by incorporating 
the total isospin T=l and T=2 states in the partial wave basis. This operation allows to 
account both for the effects due to the charge independence and charge symmetry breaking 
in the strong part of the nuclear interaction. 



III. RESULTS 

In this study proton scattering on 3 H nuclei for incident (laboratory) proton energies up 
to E p =1 MeV, that is below n- 3 He threshold, is considered. This energy region is very 
delicate due to the presence of the first excitation of the a particle, physically situated, 
just above p- 3 H threshold, see Fig. [U Subthreshold scattering cross section is very sensitive 
to the precise position of this resonance. Actually, the width of the resonance is strongly 
correlated with its position relative to the p- 3 H threshold. If this state is slightly overbound 
the resonance peak in excitation curve is naturally shifted to lower energies and becomes 
more narrow. In case of underbinding one will have much broader resonance, reflected in 
the flat excitation curve j^(E,6) as a function of energy. Very accurate description of the 
resonant scattering is therefore required. In particular, special care should be taken of charge 
symmetry breaking terms in order to properly separate the p- 3 H and n- 3 He thresholds. Our 
scattering calculations, which do not restrain the total isospin T, fully accounts for the 
isospin symmetry breaking and thus are perfectly suited. 

Several very different realistic nuclear Hamiltonians have been tested in this work. Those 
included local configuration space potential of Argonne group AV18 [19], non-local configu- 



ration space potentials INOY 20j and ISUJ 2l|, as well as chiral effective field theory based 
potential of Idaho group derived up to next-to-next-to-next-to-leading order (I-N 3 LO) [22 ]. 
Urbana three-nucleon interaction was also used in conjunction with AV18 and I-N 3 LO po- 
tentials. The parametrization, commonly called UIX [23], of this three-nucleon interaction 
was used together with AV18 model for NN interaction; on the other hand the two-pion ex- 
change parameter has been assigned slightly different value A27r=-0. 03827 MeV when using 
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this force in conjunction with I-N 3 LO NN potential. The last adjustment allowed to fit the 
triton binding energy to its experimental value in a model we refer I-N 3 LO+UIX* in that 
follows. 

Partial wave expansion is one of the most common ingredients in few-nucleon calculations. 
This expansion converges rather fast for a low energy system, since the effective centrifugal 
terms grow rapidly with the angular momentum. Nevertheless, due to multitude of the 
degrees of freedom it represents, the partial wave basis for four nucleon system becomes 
considerable and requires thousands of amplitudes to achieve numerically converged results. 
Therefore partial wave convergence turns to be a central issue for any calculation claiming 
numerical accuracy. In this study numerical convergence was studied as a function of j max = 
max(j x ,j y ,j z ) for amplitudes K and j max = max(j x ,j y ,l z ) for amplitudes H, where j x , j y , 
j z and l z are the partial angular momenta as explained in Q] . In Table [J convergence for 
p- 3 H scattering length calculations with INOY potential are presented. As one can see 
convergence pattern is not regular, however final results seem to be converged better than 
at 0.5% level. 

TABLE I: Convergence of p- 3 H scattering length calculations for INOY potential. Scattering 
lengths are presented as a function of the maximal value of the partial angular momenta (j ma x) 
allowed in the calculations. 



jmax i?( 3 H) ao ai 

1 8.049 1.070 5.051 

2 8.402 -30.65 5.583 

3 8.481 -37.01 5.452 

4 8.482 -36.87 5.377 

5 8.483 -37.16 5.371 

6 8.483 -37.35 



In the first three columns of the Table [TT] calculated ground state binding energies for 
three-nucleons and alpha particle are given. These values perfectly agree with the ones 
obtained using other ab-initio methods. The obtained results for triplet and singlet p- 3 H 
scattering lengths are also provided in Table [Til These values are estimated with some small 
error, which is due to the fact that numerical calculations have been performed for incident 
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proton energies E p > 7 keV and then extrapolated to get corresponding scattering lengths 
at E p = 0. One can see important spread of the model predicted values, in particular for 



TABLE II: Different nuclear model predictions for bound state energies of triton, 3 He and 4 He in 
MeV together with p- 3 H scattering lengths in fm. 



Model 


B( 3 H) 


B( 3 He) 


B( 4 He) 


a (p - 3 H) 


ai(p- 3 H) 

_L \± / 


ISUJ 


8.482 


7.718 


28.91 


-35.5(2) 


5.39(1) 


INOY 


8.483 


7.720 


29.08 


-37.4(2) 


5.37(1) 


AV18 


7.623 


6.925 


24.23 


-15.5(1) 


5.79(1) 


AV18+UIX 


8.483 


7.753 


28.47 


-23.6(2) 


5.47(1) 


I-N 3 LO 


7.852 


7.159 


25.36 


-19.6(2) 


5.85(1) 


I-N 3 0+UIX* 


8.482 


7.737 


28.12 


-22.9(2) 


5.59(1) 


Exp. 


8.482 


7.718 


28.30 







the resonant singlet scattering length. As already explained this value is determined by the 
position of the a particle excitation relative to the p — 3 H threshold, which are separated 
roughly by ^ h 2 /(2fial) <100 keV (here /i is reduced mass of the p- 3 H system). Therefore 
only a few keV variation of the resonance position relative to p- 3 H threshold can shift singlet 
scattering length by as much as 1 fm. Though it is difficult to establish correlation laws 
for attractive four nucleon system, one should stress sensitivity of the scattering length to 
the prediction of the ground state binding energy of the alpha particle. As one can see the 
170 keV reduction in alpha particle binding energy obtained when passing from INOY to 
ISUJ models, permits to reduce the resonant scattering length by almost 2 fm. This is not 
surprising however, once one shifts alpha particle ground state relative to p- 3 H threshold - 
its excitation is also shifted in the same direction. 

The p- 3 H triplet channel (J 71 " = 1 + ) is repulsive, while variation of the p- 3 H triplet scat- 
tering length is rather moderate: this scattering length tends to decrease as three-nucleon 
binding energy increases. Correlation is not perfect however, in particulary results for mo- 
mentum space I-N 3 LO potential tends to break the correlation. Probably it is due to the 
fact that for this potential we have used point-Coulomb repulsion, while the configuration 
space potentials employ screened Coulomb interaction having the same parametrization as 
the one encoded in EM part of the AV18 proton-proton potential. 
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In the Fig. [2] calculated p- 3 H differential cross section at 8 cm =l20° as a function of proton 
laboratory energy is compared to the experimental data. At very low energies, below 100 
keV, scattering cross section diverges due to the effective Coulomb repulsion between the 
proton and triton nucleus. This behavior, which is easily described analytically, hides com- 
pletely strong interaction effects. Only beyond 100-150 keV nuclear scattering amplitudes 
are pronounced, revealing the clear resonant behavior of the singlet channel (J 77 = + ), even 
though this wave is suppressed by the statistical factor 3 compared to the triplet (J 71- = 1 + ) 
one. The resonant peak is in particular neat for INOY and ISUJ models, which places 
it too close to the p- 3 H threshold together strongly underestimating its width. The res- 
onance width naturally increases as resonance moves further into continuum - therefore 
other model predicted excitation functions (j^(E, 9 = 120°)) are much natter than INOY or 
ISUJ ones. Probably the best description of the experimental data in the resonance region 
(E p ^200-600 keV) is achieved by I-N3LO model. The AV18 Hamiltonian slightly underpre- 
dicts experimental scattering length, its excitation function is too flat and elastic scattering 
cross section is underestimated below 600 keV. Inclusion of Urbana type three nucleon force 
makes I-N3LO+UIX* and AV18+UIX model results almost indiscernible, while the elastic 



S 



cross section is visibly overestimated in the resonance region. Th e ag reement between two 



model predictions is probably coincidental, as Hale and Hofmann|12| demonstrated strong 
sensitivity of the p- 3 H cross section to the 3NF parameters. Qualitatively the same behav- 
ior is observed, when analyzing the differential cross section as a function of the scattering 
angle. It is demonstrated in the left pane of the figure [3] for the incident protons of 300 keV. 




FIG. 3: Calculated p- 3 H differential cross sections are compared with experimental data of |28| 
for incident proton energies of 0.4 MeV (left pane) and 0.9 MeV (right pane). If at 0.4 MeV 
contribution of P-waves is negligible, they are necessary at 0.9 MeV in order to reproduce bending 
in the experimental curve. On the other hand if S-wave contributions to the elastic scattering cross 
section differ from model to model at 0.3 MeV, they merge the same curve at 0.9 MeV. 

Beyond the resonance region the S-wave cross sections for all models coincide. The AV18 
prediction integrates into the joint curve the last, only for incident proton energies around 
1 MeV. This is due to the fact that AV18 places the resonance too far into continuum, 
overestimating its width and resulting extended in energy resonance region. The 4 He spectra 
contains several P-wave resonances just above n- 3 He threshold, see Fig.[TJ These resonances, 
and in particular narrow J n = 0~ one, extend also into the region below n- 3 He threshold. 
Therefore description of the scattering cross section already beyond E p ?a400 keV requires 
inclusion of the negative parity states. This makes calculations largely involved, thats why 
calculations including P-waves have been performed only for INOY potential. These waves 
seems to improve agreement with experimental data measured close to n- 3 He threshold. 
This region is best to be studied analyzing differential cross section as a function of the 
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scattering angle: it is done in the right pane of the figure [3] for incident protons of 900 
keV, which is only ~100 keV below n- 3 He threshold. One can see that at this energy S-wave 
contribution to the differential cross section is almost the same for all the models considered. 
As demonstrated for INOY model: P-waves, having in particular strong contribution from 
the J 71 " = 0~ channel, are required in order to reproduce the bending of the experimental 
cross section curve at the backward angles. Forward scattering cross section is however 
slightly underestimated by this model, indicating that negative parity resonances are slightly 
displaced to higher energies by this interaction. The similar observation have been made by 



Deltuva et al. 



241 ] for low energy n- 3 He scattering. It is obvious that NN P-waves should have 



a strong impact for the negative parity states of the nuclear system and therefore stronger 
NN P-waves than those given by the INOY model will be favorised by the experimental 
data. 

Before closing this section possibility of using 4 He continuum as a laboratory to fine tune 
the two and three nucleon interactions should be emphasized. Indeed, the p- 3 H scattering 
at very low energies, having particular sensitivity to charge symmetry breaking terms in NN 
S-waves, is important to understand the charge symmetry breaking of the nuclear interac- 
tion. At higher energies close to n- 3 He threshold contribution of negative parity states to the 
scattering process can be well separated from the contribution of the positive parity ones. 
Negative parity states turns to demonstrate stronger sensitivity to NN P-waves. Furthermore 
if the study of the 4 He continuum is undertaken in parallel with resonant scattering in n- 3 H 
and p- 3 He systems the new frontiers opens to understand isospin symmetry of NN P-waves 
as well as isospin structure of the three-nucleon interaction. In view of such perspectives im- 
portance of low energy experiments in four nucleon system should be stressed. In particular 
the nucleon scattering on tritium data is very scarce, while the tritium experiments have 
been completely abandoned for the last 30 years. The nice exception is recently performed 
proton-tritium differential cross section measurement in Fudan University (Shanghai) 29] 
for proton incident energies from 1.4 MeV to 3.4 MeV and detectors fixed at the backward 
laboratory angle of 165°. Hopefully this activity is continued for lower energy protons as 
well as extended to use the neutrons as projectile. 
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IV. CONCLUSION 



In this work the proton elastic scattering on tritons is studied for incident proton ener- 
gies below 1 MeV using Faddeev-Yakubovski equations in configuration space and by fully 
accounting for isospin breaking effects in nuclear interaction. Converged results are ob- 
tained for six qualitatively different realistic nuclear Hamiltonians. However none of the 
tested Hamiltonians have been able to reproduce the shape of the nearthreshold J n = 0+ 
resonance together with the ground state binding energy of the alpha particle. 

It has been demonstrated that 4 He continuum is an ideal laboratory to fine tune the 
two and three nucleon interaction models and the particular role this system can play in 
understanding the charge symmetry breaking. Recent advances in computational techniques 
allow very precise calculations for four-nucleon scattering problem at low energies, therefore 
new four-nucleon scattering experiments and in particular ones concerning nucleon scattering 
on tritium are strongly anticipated. 
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